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Preface 


The  purpose  of  this  thesis  was  to  analyze  the  causes  of 
subharmonic  aliasing,  develop  an  analytic  model  of  this 
effect, and  determine  its  characteristics  on  the  AFTI/F-16 
Digital  Flight  Control  System  (DFCS). 

Many  hybrid  (SIMSTAR)  computer  simulations  were 
performed  to  verify  that  subharmonic  aliasing  existed  and 
determine  its  effects  on  the  asynchronous,  triply  redundant 
digital  flight  control  system.  After  weeks  of  trying  to 
model  the  effect  in  the  time  domain,  I  was  steered  toward 
the  Fourier  transform  and  the  frequency  domain  by  Lt  Col 
Lewantowlcz.  Within  days  a  model  was  formed. 

In  performing  the  simulations  and  modelling  the 
phenomenon,  I  received  a  great  deal  of  help  from  others.  I 
would  like  to  thank  my  thesis  advisor,  Lt  Col  Lewantowlcz, 
for  his  guidance  in  helping  me  model  subharmonic  aliasing. 

I  would  also  like  to  thank  Bob  Ewing  for  his  help  with 
SIMSTAR  and  for  providing  a  work  area  used  for  the  many 
hours  of  programming  necessary  to  complete  this  research, 
and  Don  Smith  (EAI)  for  his  persistance  in  keeping  the 
SIMSTAR  "up".  Finally,  I  would  like  to  thank  my  family, 
Jeremy,  Nathan,  Bethany,  and  especially  my  wife  Cindy,  for 
their  patience,  love,  and  understanding  throughout  the 
entire  eighteen  months  at  AFIT. 


David  M.  Thomas 
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Abstract 


The  purpose  of  this  research  is  threefold.  First, 
determine  the  cause  of  subharmonic  aliasing,  described  by 
the  AFTI/F-16  engineers  as  "the  creation  of  uncorrelated 
low  frequencies  whenever  a  subharmonic  of  the  sample 
frequency  is  input  into  the  system".  Second,  model  the 
subharmonic  aliasing  effect,  so  that,  by  knowing  only  the 
input  frequency  and  the  system  sample  rate  the  output 
characteristics  can  be  calculated.  And  third,  demonstrate 
by  simulation  the  effect  of  input  and  output  filters  on  the 
subharmonic  alias,  and  the  effect  of  signals  in  the 
subharmonic  range  (Wg/lO  <  Wq  <  W5/2)  on  the 
interchannel  difference  and  the  software  rate  limiter. 

The  model  determined  that  subharmonic  aliasing  is  the 
result  of  impostor  frequencies  (much  like  aliasing)  being 
introduced  into  the  output  signal  by  the  sampling  process. 
Well  defined  subharmonic  aliases  occur  due  to:  impostor 
frequencies  and  a  phenomena  known  as  apparent  low  frequency 
surge,  which  occurs  when  the  input  frequency  is  nearly  an 
integer  multiple  (>1)  of  the  impostor  frequency. 

Simulations  demonstrate  the  effectiveness  of  input 
analog  low  pass  filters  for  attenuating  signals  in  this 
frequency  range.  Unfortunately,  according  to  simulations. 


Interchannel  differences  may  exceed  Input  or  Output  selector 
monitor  threshold  levels  for  frequencies  in  the  subharmonic 
range.  Simulations  also  show  that  it  is  the  combination  of 
well  defined  subharmonic  aliases,  the  high  frequencies 
(relative  to  the  sample  rate)  of  the  input  signal,  and  the 
nonlinear  characteristics  of  the  software  rate  limiter  which 
form  the  low  frequency  output  from  the  rate  limiter. 

In  conclusion,  since  subharmonic  aliasing  occurs  for 


frequencies  which  are  high  relative  to  the  sample  rate, 
careful  selection  of  the  sample  rate  (and  computational 
frame  rate)  should  reduce  this  undeslreable  effect. 
Additionally,  since  the  phenomenon  is  actually  the 
combination  of  two  high  frequencies,  the  use  of  analog 
lowpass  filters  at  the  input  does  attenuate  the  subharmonic 
aliasing  effect. 


SUBHARHONIC  ALIASING  AND  ITS  EFFECTS  ON 


THE  AFTI/F-16  DIGITAL  FLIGHT  CONTROL  SYSTEM 


I .  Introduction 

Vith  the  advent  of  smaller,  faster  computers  came  the 
pioneering  of  digital  flight  control  systems  (DFCS)  to  steer 
future  aircraft  (1:1).  The  Advanced  Fighter  Technology 
Integration  (AFTI)  F-16  is  such  an  aircraft.  Its 
asynchronous,  triple-redundant,  digital  flight  control 
system  was  designed  to  meet  all  flight  control  needs 
previously  performed  by  a  quad-redundant  analog  system. 
Because  the  AFTI/F-16  is  completely  fly-by-wlre  control 
(there  are  no  mechanical  backups  to  the  flight  control 
system),  considerable  effort  was  made  to  ensure  reliability, 
safety,  and  performance.  In  case  of  hardware  failure  within 
the  DFCS,  an  advanced  strategy  for  redundancy  management  was 
developed  to  both  detect  and  isolate  DFCS  failures. 

Redundant  elements  (sensors, computers,  etc)  were 
designed  into  the  system,  and  to  monitor  the  inputs  and 
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outputs  o£  these  elements  an  advanced  strategy  for 
redundancy  management  was  developed,  in  such  a  design, 
comparison  monitoring  and  voting  is  used  to  detect,  isolate, 
and  remove  failures.  These  redundant  elements  and  the 
strategy  for  redundancy  management  will  be  described  briefly 
in  Chapter  2. 

Two  problems  are  Inherent  to  the  AFTI/F-16  DFCS: 
interchannel  differences  caused  by  time  skew  and  subharmonic 
aliasing.  The  first  problem  is  unique  to  the  AFTI/F-16 
since  its  DFCS  is  asynchronous.  In  effect,  the  asynchronous 
operation  causes  Interchannel  differences,  because  of  time 
skew,  between  the  three  redundant  channels.  The  inter¬ 
channel  difference  may  be  Interpreted  by  the  fault  detection 
and  isolation  algorithms  as  a  failure.  This  problem  is 
discussed  extensively  in  Lt  Gursel  Serpen's  thesis  "Failure 
Detection  and  isolation  for  an  Asynchronous  Digital  Flight 
Control  System"  (2). 

Problem 

The  second  problem,  the  main  topic  of  my  thesis,  is 
subharmonic  aliasing.  Subharmonic  aliasing  is  a  phenomenon 
described  by  the  AFTI/F-16  engineers  as  the  creation  of 
uncorrelated  low  frequencies  whenever  a  subharmonic  of  the 
sampling  frequency  (fg)  is  input  into  the  system  (1:10). 

This  study  shows  that  it  is  actually  a  phenomena  which 


effects  all  sample  data  systems  that  do  not  have  Ideal  anti¬ 
aliasing  filters.  Subharmonic  aliasing  is  most  prominent 
for  sinusiodal  Inputs,  whose  frequencies  f^  (Hz)  are  less 
than  t^/2  (the  Nyquist  rate)  but  greater  than  fg/10. 

Figures  l(a  &  b)  show  two  examples  of  subharmonic  aliases. 

The  input  is  a  sine  wave,  the  output  is  the  input  sampled  at 
64  samples/second  (the  OFCS  sample  rate).  A  main  concern  of 
this  thesis  is  that  the  AFTI/F-16  has  a  number  of  structral 
modes  within  a  range  of  frequencies  labeled  the  subharmonic 
alias  range,  fs/10  Hz  <  f^  <  fg/2  Hz,  which  for  fs=  64  samples/ 
second  is  6.4  Hz  <  f^  <  32  Hz  (see  Table  I). 

Because  of  the  asynchronous  DFCS  architecture  of 
the  AFTI/F-16,  the  high  frequency  component  of  the 
subharmonic  aliased  signal  may  cause  large  interchannel 
differences  between  redundant  channels  and  may  appear  as  a 
disagreement  to  the  signal  monitors.  This  disagreement  may 
trigger  the  fault  detection  and  isolation  algorithm,  which 
in  turn  may  vote  a  sensor  as  failed  when  no  failure  actually 
exists.  Abnormally  high  Input  and  Output  selector  monitor 
threshold  levels  are  used  to  tolerate  this  in  the  AFTI/F-16. 
(This  paragraph  directly  contradicts  the  AFTI/F-16  engineers 
assessment  that  it  is  the  low  frequency  component  that  would 
cause  the  Interchannel  differences.  The  problem  is 
explained  in  detail  in  Chapter  5.) 
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Table  I.  AFTI/F-16  Structural  Modes 


Symmetr ic 


Frequency  (Hz) 


4.59 


Antlsymmetr ic 


11.26 


19.235 


20.463 


21.186 


35.592 


37.293 


5.798 


8,343 


12.753 


15.887 


Damping  (g) 


.023 


.040 


.042 


.056 


.040 


.056 


.020 


.023 


.022 


.044 


.060 


19 . 177 


22.008 


.046 
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The  major  assumption  for  this  thesis  is  that,  although 
the  control  laws  and  aircraft  dynamic  model  are  not 
simulated,  the  sample  rates  and  inputs  are  simulated  and 
this,  adequately  demonstrates  subharmonic  aliasing  and  its 
effects  on  the  DPCS. 


Scope 


This  thesis  describes  in  detail,  mathematically  and  by 
simulation,  subharmonic  aliasing.  It  shows  how  frequencies 
in  the  "subharmonic  aliasing  range"  effect  the  interchannel 
differences  between  the  Digital  Plight  Control  Computers 


because  of  time  skew  and  how  the  actual  threshold  levels  for 
the  Input  S/M  and  the  Output  S/M  can  be  violated.  Plnally, 
effects  of  the  rate  limiter  are  demonstrated  and  two  config¬ 
urations  of  the  rate  limiter  placement  within  the  DPCS  are 
cons idered . 

In  all  hybrid  computer  (SIHSTAR)  simulations.  Inputs 
are  generated  by  an  Internal  harmonic  oscillator  (sine  wave 
generator)  rather  than  external  inputs  to  limit  noise 
effects  and  for  ease  of  parameter  control. 


The  following  approach  la  used  to  Investigate  the 
"Effects  of  Subharmonic  Aliasing  on  the  AFTI  F-16  Digital 
Flight  Control  System". 

1.  A  simple  sampling  system  Is  simulated  on  the 
hybrid  computer  (SIMSTAR)  to  recreate  the  subharmonic 
aliasing  effect  as  described  by  the  AFTI/P-16  engineers. 

2.  A  mathematical  model  Is  developed  and  used  to 
predict  the  magnitude  and  frequency  of  the  subharmonic 
aliased  signal.  A  digital  computer  simulation  package 
Is  used  to  verify  the  model. 

3.  The  maximum  Interchannel  difference  is 
mathroatlcally  computed  and  graphed  for  the  frequencies  in 
the  subharmonic  range.  Hybrid  and  digital  computer 
simulations  are  used  to  verify  the  computations. 

4.  A  software  rate  limiter  is  simulated  using 

a  FORTRAN  subroutine  In  a  hybrid  computer  to  determine  Its 
effect  on  subhamonlc  aliased  signals. 

5.  Means  to  reduce  subhamonlc  aliasing  are  examined 

6.  Follow-on  work  Is  suggested. 

The  Resources  required  to  accomplish  this  thesis  are 
listed  In  Appendix  C. 


II .  BACKGROUND 


Xlie  AFT?/Frl6  Ptattat  Fiiq.ht  goptgol  gyst^pt 

The  APTI/P-16  Digital  Flight  Control  System  (DFCS)  is  a 
full  authority,  triplex,  digital  fly-by-wire  flight  control 
system  that  Includes  provisions  for  stringent  reliability 
and  fail-safe  operational  requirements  (3:4-1).  The  triplex 
configuration  consists  of  three  identical  flight  control 
computers  (FLCC)  which  operate  asynchronously  with  respect 
to  each  other.  They  are  designed  to  operate  mutually 
independent  of  each  other  for  increased  system  reliability. 
Each  FLCC  treats  the  other  two  as  remote  sensors  rather  than 
as  an  external  computational  source,  and  each  FLCC  forms  an 
independent  view  as  to  the  state  of  the  system.  Therefore, 
it  is  possible  for  one  computer  to  differ  from  the  others 
with  respect  to  what  is  failed  in  the  system  or  in  what  mode 
the  system  is  operating  (3:1-3). 

The  overall  DFCS  is  Illustrated  in  Figure  2.  The  three 
identical  FLCC's  provide  computational  power  for  control  law 
computation  and  redundancy  management.  A  fourth  unit  (not 
shown),  the  Actuator  Interface  Unit  (AIU),  houses  dual- 
redundant  analog  functions  associated  with  the  leading-edge 
flap  actuator  drives  (3:4-23). 


The  Computer 


The  highlights  of  the  computer  architecture  are 
presented  here,  while  detailed  Information  Is  contained  In 
Reference  3.  Each  FLCC  contains  the  following  eight  major 
functional  sections: 

1.  Analog/discrete  input  signal  conditioning  and 
multiplexing 

2.  Analog/discrete  output  signal  conditioning  and 
demultiplexing 

3.  Digital  computation  and  data  conversions 

4.  Intercomputer  data  link  Interface 

5.  Multiplex  bus  (MIL  -  STD  -  1553)  interface 

6.  Dedicated  failure  logic 

7.  ISA  drives,  associated  switching,  and  warning 
light  drives 

8.  Analog  Independent  backup  unit  (IBU) 

The  digital  portion  of  the  computer  architecture  is 
composed  of  the  following  five  basic  sections: 

1.  Central  Processing  Unit  (CPU) 

2.  Memory  (CPU  Random  Access  Memory) 

3.  Input/output  controller  (IOC) 

4.  Intercomputer  data  link  receiver  controller 

5.  The  1553  Multiplex  bus  Interface  controller 

In  this  architecture,  the  processor  communicates  only  with 
the  memory  and  performs  no  Input/output  (I/O)  operations 


except  for  those  critical  discretes  that  are  mapped  into 
memory  (3:4-29). 


Input/Output  Controller 

The  I/O  operations  ate  controlled  by  the  Input/Output 
Controller  (IOC)  which  minimizes  the  CPU  real  time  pro¬ 
cessing  requirements.  The  IOC  supervises  and  controls  the 
operation  of  the  Input/Output  Multiplexer/Demultiplexer 
system,  the  analog  to  digital  (A/D)  and  digital  to  analog 
(D/A)  converters,  the  Intercomputer  Data  Link  Transmitter 
and  the  Input/Output  Discrete  Registers.  The  IOC  is  micro- 
programmable  to  allow  a  random  sequence  of  I/O  operations 
(up  to  1000  unique  operations),  and  in  the  autonomous  mode 
its  operation  is  completely  transparent  to  the  processor. 
For  normal  (non-failure)  conditions,  the  processor  retains 
supervisory  authority  over  the  operation  in  that  it  can 
vector  the  controller  to  specific  prepro-grammed  I/O 
sequences,  and  it  receives  an  indication  from  the  IOC  that 
the  requested  sequence  has  been  completed.  This  provides 
reasonable  control  over  the  latency  of  data  to  be  used  for 
the  control  law  or  redundancy  management  functions.  The 
individual  operations  of  the  command  sequence,  however, 
remain  transparent  to  the  processor.  The  IOC  selects  the 
I/O  device.  Initializes  conversions  and  transfers  data  to/ 
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from  memory  according  to  a  preprogrammed  file,  hence  the 
processor  needs  only  to  access  the  memory  for  I/O  opera¬ 
tions.  For  detected  failures  of  the  channel  (e.g.,  the 
processor),  the  dedicated  failure  logic  overrides  the  pro¬ 
cessor  control  and  vectors  the  I/O  controller  to  a  special 
file  which  maintains  the  unit  in  a  free-running  mode  until 
the  failure  condition  is  removed.  This  free  running  auton¬ 
omous  mode,  allows  the  cross-channel  communication  of  sensor 
data  (and  fault  status)  to  be  continued  to  the  other  two 
channels  of  the  triplex  system  even  after  a  processor 
failure  (4:1429). 


Monitoring.  Voting,  and  Failure  Management 

Most  of  the  monitoring,  voting,  and  redundancy  manage- 
gement  functions  within  the  flight  control  computer  are 
implemented  in  software.  Voter/monitor  planes  are  provided 
in  each  channel  for  redundant  sensor  inputs  by  use  of  digit¬ 
ally  cross-strapped  data  from  the  other  two  channels. 
Cross-channel  comparison  monitoring  of  sensor  and  computed 
data  is  the  primary  method  of  fault  detection  with  in-line 
self-testing  used  only  to  Isolate  a  failed  channel  following 
a  second  similar  failure  detected  by  the  cross-channel 
monitoring.  This  method  achieves  a  probability  of  loss-of- 
control  (PLOC)  of  KT^  by  the  use  of  only  three  channels  of 


computation  (5:1222). 


Provisions  are  included  to  assure  that  a  failed  branch 
is  prohibited  from  contributing  to  the  control  of  the  air¬ 
craft.  Since  a  malfunctioning  computer  cannot  be  expected 
to  recognize  its  own  errors,  special  dedicated  hardware 
islncluded  to  permit  the  unfailed  processors  to  disconnect 
the  malfunctioning  computers  outputs.  For  a  second  like- 
failure  control,  each  remaining  processor  can  alarm  the 
failure  and  require  a  response.  The  AFTI/F-16  uses  this 
approach  since  Internal  self-testing  is  relied  on  for  action 
in  resolving  second  failures  rather  than  empower  a  single 
branch  to  disconnect  another  branch.  If  the  second  like- 
failure  cannot  be  resolved  by  self-testing,  automatic 
reversion  to  an  analog  backup  control  system  is  provided 
(3:4-41) . 

Input  Selector  Monitor  Plane  ( Input  S/M) .  The  function 
of  the  input  selector  monitor  is  to  detect  failures  of 
analog  Inputs  from  the  three  redundant  FLCC’s  and  determine 
a  good  analog  input  for  use  by  the  control  law  function. 
Failure  detection  is  accomplished  by  use  of  cross-channel 
monitoring  techniques  (6:3-9). 

Note : 

Self  -  refers  to  the  FLCC  performing  the  calculation. 

Left  and  Right  -  refer  to  the  "other  two"  FLCC's. 

Analog  inputs  are  monitored  in  the  following  manner. 

The  current  value  of  a  sensor  input  is  read  from  the  self 
IOC  random  access  memory  (IOC  RAM)  (Figure  3)  and  the  latest 


value  from  the  left  and  right  serial  data  link  receiver 
(SDLR)  RAM  areas  respectively.  If  the  three  values  differ 
by  no  more  than  the  allowable  interbranch  trip  level  for 
this  parameter  (15%  of  full  scale  for  most  parameters),  then 
the  average  of  these  three  values  is  calculated  and  stored 
in  the  CPU  RAM  area  for  later  utilization  by  the  control  law 
functions.  If  one  of  the  three  input  values  differs  from 
the  other  two  an  amount  greater  than  the  trip  level,  the  out 
of  tolerance  signal  is  rejected  by  the  monitor  and  the 
selector  algorithm  calculates  the  average  of  the  remaining 
two  good  values.  This  type  of  algorithm  is  known  as  a  good 
channel  averager  (3:1223). 

If  this  out  of  tolerance  condition  persists  for  more 
than  seven  major  computer  frames  (each  frame  is  15.625  msec) 
the  failure  is  considered  a  hard  failure  and  a  call  is  made 
to  the  Failure  Management  function  to  record  the  failure  and 
to  perform  further  failure  analysis. 

Output  Selector  Monitor .  The  function  of  the  Output 
Selector  Monitor  (Output  S/M)  is  to  monitor  the  Total 
Computed  Output  (TCO)  (Figure  4)  or  surface  commands  to 
detect  computational  failures  and  to  ensure  that  the  three 
coil  current  drives  to  each  ISA  servo  valve  are  always 
within  the  ISA  hydraulic  trip  level  (10%  to  30%  of  full 
scale  coil  current)  so  long  as  the  command  path  is  failure 
free.  The  design  of  the  output  monitor  is  similar  to  the 
Input  S/M  design  and  Includes  three  persistence  counters, 
one  totalizer  counter,  separate  first  fail  and  second  fail 


maximum  persistence  values,  and  a  fixed  allowable  cross¬ 
channel  difference.  In  order  to  minimize  latency  through 
the  selector.  It  Is  necessary  that  the  selector  pick  the 
output  of  the  computer  that  Is  actually  driving  the  surface. 
When  the  system  Is  operating  from  the  primary  hydraulic 
system  (servo  valve  3),  minimum  latency  Is  achieved  by 
selecting  FLCC  C  to  drive  the  system  (6:5-39).  (For 
configuration  management  purposes  the  FLCC's  are  Identified 
as  FLCC  A,  FLCC  B,  and  FLCC  C). 

Failure  Management .  The  Failure  Management  function 
provides  centralized  management  of  the  DFCS  failures.  This 
function  checks  the  status  of  the  DFCS  by  referring  to  the 
Device  Status  Table  (DST),  it  records/updates  the  device 
failure  on  the  DST  to  remove  the  failed  element  from  the 
redundancy  scheme,  and  it  maintains  a  count  of  low  level 
device  failures  to  determine  if  a  high  level  device  is 
actually  failed  (6:11-1,11-19). 


Data  Acquisition 

Several  key  factors  are  Important  in  developing  the 
overall  data  acquisition  scheme  for  the  AFTI/F-16  DFCS. 
Only  one  of  these  factors  Is  presented  here.  The  major 
factor  that  determines  the  rate  at  which  data  must  be 


acquired  is  the  rate  of  change  of  the  input  signal  with 


respect  to  time.  In  general,  if  the  rate  of  change  of  the 
input  is  high  relative  to  the  sampling  rate,  a  large  inter¬ 
channel  difference  is  generated  in  otherwise  identical  sig¬ 
nals  when  the  data  is  exchanged  between  computers  (6:2-1). 

For  example,  the  largest  pilot  stick  pitch  input  (s) 
expected  is  0  to  9  (full  scale)  in  0.1  seconds,  in  other 

s 

words  an  input  rate  of  90  g/second  (s).  The  maximum  sample 
period  (T)  allowed  to  keep  the  interchannel  difference  (e) 
less  than  5%  can  be  calculated  by  the  simple  equation 

T  =  e^3x/s  (1) 

=  0.45  g/(90  g/sec) 

e  0.005  second 

where 

T  =  period 

e^gjj  =  maximum  interchannel  difference 
=  9  g  *  0.05 
=  0.45  g 

s  =  maximum  signal  rate  of  change  (time  derivative) 

=  9  g  /  O.lsec 
=  90  g/sec 

Since  the  DFCS  has  a  computational  frame  period  of  1/64 
seconds  or  15.625  msec  (6:2-8),  for  ease  of  implementation, 
a  multiple  of  this  period  is  required  for  input  data 
acquisition.  The  input  period  selected  is  3.90625  msec 


(15.625/4)  or  256  samples  per  second. 

As  shown  in  Figure  5  (6:2-3)  this  rate  is  adequate  even 
when  combined  with  anti-aliasing  filters  or  prefilters  with 
cutoff  frequencies  at  u  =  100,  50  and  10  rad/sec.  However, 
as  break  frequencies  decrease,  excessive  phase  lags  are 
Introduced  resulting  in  reduced  gain  and  phase  margins  and 
poor  handling  characteristics  (6:2-2). 

Since  interchannel  differences  are  not  presently 
reduced  by  voter  algorithms  (as  are  sensor  bias  and  calibra¬ 
tion  gain  errors),  minimizing  these  errors  with  an  adequate 
sampling  rate  is  critical  for  the  AFTI/F-16  (6:2-2).  (See 
Ref  2,  Lt  Gursel  Serpen's  thesis,  for  a  method  to  reduce  the 
interchannel  dlfferece  by  estimating  the  time  skew.) 

Rate  Limiters 

Software  rate  limiters  are  mechanized  in  the  surface 
command  mixer  just  after  the  Output  S/M  plane  (Figure  4)  for 
each  control  surface  except  for  the  leading  edge  flap.  The 
rate  limits  are  set  to  limit  the  signal  change  per  Iteration 
to  each  surface  to  no  more  than  5%  of  full  scale  deflection. 
The  purpose  of  this  rate  limiting  is  twofold.  First,  it 
acts  to  minimize  interchannel  differences  caused  by  time 
skew,  and  second,  it  acts  to  prevent  erroneous  ISA  output 
electronic  trips  due  to  large  f rame-to-f rame  command 
deviations  (6:13-23). 


Aliasing 


Aliasing  Is  a  phenomena  well  described  In  numerous 

texts  (8,9,10)  on  digital  control  and  communications  and  can 

be  found  under  headings  of  either  "aliasing"  or  "folding". 

It  occurs  when  Shannons  Sampling  Theorem  Is  violated. 

Shannons  Sampling  Theorem  simply  states: 

Let  x(t)  be  a  bandllmlted  signal  with  X(u)  =  0 
for  IwK  Wjp.  Then  x(t)  Is  uniqully  determined  by 
its  samples  x(nT),  n=0,  +1,  +2,...  if 

Ws  >  2w^  (2) 

where 


Ug  -  2ic/T  sample  radian  frequency 
T  s  sample  period 

=  maximum  radian  frequency 

Given  these  samples,  x(t)  can  be  reconstructed 
by  generating  a  periodic  impulse  train  in  which 
successive  impulses  have  amplitudes  that  are 
successive  sample  values.  This  impulse  train  Is 
then  processed  through  an  Ideal  low-pass  filter  with 
gain  T  and  cutoff  frequency  greater  than  and 
less  than  (Wg  -  w^)  .  (7:519). 


This  theorem  Is  violated  when 


“m  > 

where  «s/2  is  refered  to  as  the  Nyqulst  rate,  the  rate  which 
must  not  be  exceeded  by  to  prevent  aliasing. 

For  example  let 


x(t)  =  COS(WQt) 


(3) 


where 


with  Fourier  transform  X(w)  as  indicated  In  Figure  6 
(9:528).  In  this  figure,  the  Impulse  Is  distinguished  at 
from  that  at  -w^  for  convenience,  as  the  discussion 
proceeds.  Consider  Xp(u),  the  spectrum  of  the  sampled 
signal  and  focus  in  particular  on  the  effect  of  a  change  in 
the  input  frequency  Wq  with  the  sampling  frequency  fixed 
Figure  6  (b)-(e)  illustrates  several  values  for  Uq.  Also, 
indicated  by  the  dashed  line,  is  the  ideal  passband  of  a 
low-pass  filter  with  »  Ws/2.  Mote  that  no  aliasing 
occurs  in  (b)  or  (c),  since  Uq  <  Wg/2,  whereas  aliasing 
does  occur  in  (d)  and  (e).  For  each  of  the  four  cases,  th 
ideal  low-pass  filter  output  of  X|.(t)  is  given  by: 


b) 

“o 

=  «s/6; 

x^  (t) 

=  cos  («ot) 

= 

x(t) 

c) 

Wo 

«  2Wg/6; 

X 

cf 

=  COS(WQt) 

x(t) 

d) 

«o 

=  4Ws/6; 

X^  (t  ) 

=  COS(«s-Wo)t 

x(t) 

e ) 

Wo 

*  5Wg/6; 

X^  ( t  ) 

=  COS(W5-«o)t 

x(t) 

When  aliasing  does  occur,  the  original  input  frequency 


Wq  ta)ces  on  the  identity  or  "alias"  of  a  lower  frequency 
«.  (7:527-529) 


-  U) 


a 


Figure  6.  Effect  In  the  Frequency  Domain  of  Adequate 
Sampling  (b)  &  (c),  and  Undersampling  (d)  &(e) 
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III.  Simulation  Model 


Figure  7  is  a  model  of  the  pertinent  hardware/software 
relationships  of  the  AFTI/F-16  digital  flight  control  system 
which  are  used  to  study  the  effects  of  frequencies  in  the 
subharmonic  range.  For  this  study  the  subharmonic  range  is 
defined  as  the  frequencies  Wq  where  w^/lO  <  <  u^/2. 

Inputs  are  generated  with  a  sign  wave  generator  that  is 
internal  to  the  hybrid  (SIMSTAR)  computer.  Values  for  the 
anti-aliasing  filters  and  prefilters,  if  used  in  the 
simulation,  are  from  Table  II  which  is  a  list  of  anti¬ 
aliasing  and  pre-filters  used  by  the  AFTI/F-16  (14).  An 
input  sample  rate  of  256  samples/second  is  used  to  simulate 
the  sensor  sample  rate  of  the  aircraft.  The  second  sample 
rate  of  64  samples/  second  simulates  the  computational  frame 
rate  of  the  FLCC  described  in  Chapter  II. 

Chapters  4,  5,  and  6  implement  different  parts  of  this 
model  to  show 

1.  Vhat  subharmonic  aliasing  is. 

2.  How  frequencies  in  the  subharmonic  range  affect  the 
magnitude  of  the  interchannel  difference. 

3.  How  and  why  the  rate  limiter  'follows'  the  low 


frequency  component  of  the  subharmonic  aliased 
signal  and  where  the  rate  limiter  should  be 
(before  or  after  the  Output  Selector  Monitor  (S/M). 


Filter 


Input 

S/M 


Output 

S/M 


Rate 

Limiter 


256  samples/sec 
Sensor  Sample 
Rate 


64  aaaplea/sec 

Ccnnputatlonal 
Frame  Rate 


Figure  7 .  Hardware/Software  Model 


Table  II.  Selected  Analog  Anti-Aliasing  and  Pre-Filters 


Mode 


Stick  Input 


Normal  Accel 
(Anti -Alias ) 


Pitch  Rate 


Angle  o£  Attack 
Yaw  Rate 
Pitch  Rate 

Lateral  Accel 
(Anti-Alias ) 

Rudder  Pedal 
Throttle  Twist 
Grip 


Filter (s) 


295  *  112 

(s+295)  (s+112) 

64.2 

(Sf64.2) 


325  *  98.6 

(s+325)  (s+98.6) 


321  *  85.7 

(s+321)  (s+85.7) 

96.2 

(s+96.2) 


295  »  112 

(8+295)  (s+112) 


Impact  Pressure 
(Anti-Allas) 


96.2 

(s+96.2) 


IV, 


jLljaginq 


Subhatmonic  aliasing  is  a  phenomena  dual  to  aliasing 
and  associated  with  non-ideal  sample  data  systems,  that  is, 
systems  without  ideal  anti-aliasing  filters  with  cutoff 
frequencies  at  Wg/2,  which  use  a  hold  device  (usually  a 
zero-order  hold)  at  its  output.  The  purpose  of  this  chapter 
is  to  describe  and  model  subharmonic  aliasing.  It  also 
provides  a  number  of  examples  of  subharmonic  aliasing  and 
the  model,  provides  a  means  for  predicting  the  magnitude  and 
frequency  of  the  subharmonic  aliased  signal,  and  shows  the 
utility  of  the  model  for  prediciting  the  output 
characteristics  of  a  system. 


Description  of  Subharmonic  A1 iaslnq 

Given  the  simple  circuit  shown  in  Figure  8,  let 
x(t)  =  M*  cos(WQ*t).  The  Fourier  transform  X(w)  of  the 
input  is  sho%m  in  Figure  9(a)  and  the  frequency  spectrum  of 
the  sampled  input  X  (w)  is  shown  in  Figure  9(b).  This  shows 
that  there  are  components  of  the  input  signal  not  only  at 
but  also  at  nWg  ±  w^,  where  n  =  0,  +1,  +2,...  . 

Next,  given  the  zero-order  hold  (ZOH) 


Figure  6.  Sample-Data  System 


the  frequency  response  can  be  found  by  replacing  s  by  ju 

=  (l-exp(-Tj«)  )/jw 

which/  after  much  mathmatical  manipulation  (6)  can  be 
written  as 


®ho  *  T  *  sin ( itw/ws  )  *  exp- jx ( w/ws  )  < 

XW/Us 

The  gain  characteristics,  shown  in  Figure  9(c),  imply 
that  the  ZOH  behaves  as  a  low  pass  filter.  However,  when 
compared  to  the  gain  characteristics  of  the  ideal  filter 
(Pig.  6,  Chapter  Two),  the  amplitude  response  of  the  ZOH  is 


zero  at  w  B  Ug,  instead  o£  cutting  off  sharply  at  i^/2. 
This  provides  the  impetus  for  the  subharmonic  aliasing 
phenomena . 


The  magnitude  and  phase  characteristics  of  the  ZOH  are 
given  by  (11:134-135) 


♦  («)  =  -Ttti>  (8) 

IG^^(ju)|  =  T  IsinCzu/Wg ) I  (9) 

Ktf/ws 

Note  that  the  phase  of  the  ZOH  is  linear  with  ,  and  the 
gain  is  described  by  the  sine  function. 

When  the  composite  spectrum  of  Figure  9(a)  is  filtered 
by  the  ZOH  of  Figure  9(c),  impostor  frequencies,  shown  in 
Figure  9(d),  appear  in  the  output  signal  (12:85),  and 
create,  as  the  AFTI/F-16  engineeres  call  them,  subharmonic 
aliases . 


Examples  of  Subharmonic  Aliasing 

The  previous  section  describes  subharmonic  alaislng. 
It  is  not,  as  was  feared  by  the  AFTI  engineers,  a  problem 
to  the  AFTI/F-16  DFCS .  This  section  simply  provides  a 
number  of  examples  of  subharmonic  aliasing,  generated  by  a 
hybrid  computer  (see  Appendix  B,  SIMSTAR  PROGRAMS).  The 
simple  A/D  to  D/A  circuit  in  Figure  10  is  used  to  generate 


Figure  10.  Subharmonic  Aliasing  Simulation 

input  signals  vhich^  when  sampled,  create  impostors 
frequencies . 

For  an  input 


The  sampled  output 


r(t)  =  sin(«o*t) 

C*(s)  =  R*(5)»Z0H 
C(Z)  =  R(2)*Z0H 
c(kT)  =  sin(Wo*  k»T) 


where 


“o 

<  Ws/2 

«s 

=  the  ; 

T 

»  2i(/w 

k 

CM 

O 

H 

30 


(10) 


(11) 
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This  simulation  uses  a  sample  rate  of  64  samples/second 


(the  APTI/F-16  control  lav  computational  frame  rate). 
Figures  11-16  show  how  the  impostor  frequencies  (Ws-Wq) 
effect  the  sampled  input.  For  more  examples  on  how  the 
Impostor  frequencies  affect  the  Input  see  Appendix  A. 


The  Model 

Subharmonic  aliasing  is  modeled  using  the  information 
from  figure  9(d),  which  shows  how  an  impostor  frequency  is 
introduced  into  the  sampled  signal  due  to  non-ideal  anti¬ 
aliasing  filters  or  the  ZOH.  As  wo  becomes  small  with 
respect  to  Ws/2  (as  demonstrated  in  Figures  A1-A7,  in 
Appendix  A)  the  output  signal  approaches  a  function  of  only 
Wq  since  the  impostor  frequency  ~  (do/  which  creates  the 
subharmonic  alias,  is  attenuated  to  nearly  zero  magnitude  by 
the  filter  or  ZOH.  Given  this  information,  a  model  of  the 
sampled  signal 

c(kT)  =  M*cos(Wo*k»T)  (12) 


where 


Wq  <  Ws/2 

which  adequately  describes  the  frequency  and  magnitude 
of  the  subharmonic  alias  is  given  by 


1-1  sec-1 


Input 
12.9  Hz 


)  i  «  tr=>  t  •  i  if 


_( - H 


Output 


Figure  15. 


t  (seconds) 

Subharmonic  Aliasing  Effect  (input  12.9  Hz, 
sampled  at  64  samples/sec) 


I  nput 
9  Hz 


Output 


Figure  16 


Subharmonic  Aliasing  Effect  (input  9  Hz, 
sampled  at  64  samples/sec) 


i!!i.<!a 


c  (t)  «  M*(A*cos((iiot)-<-B*co8(  (Wg-Wo  1 

where  M  Is  the  magnitude  of  the  input  signal  prior  to 

sampling.  If  H  «  1  then  A  is  the  magnitude  of  this  signal  due  to 

attenuation  by  the  ZOH  and  B  is  the  magnitude  of  the 

impostor  frequency  due  to  the  non  ideal  filter 

characteristics  of  the  ZOH. 

Two  constraints  are  imposed  on  the  values  of  A  and  B. 

First, 

A  -t-  B  =  1  (14) 

for  all  frequencies  modeled.  This  constraint  is  imposed  so 
that  the  magnitude  M  (Eq.  13)  of  the  model  is  equal  to  the 
magnitude  M  (Eq  12)  of  the  input  signal  for  all  frequencies 
modeled.  The  second  constraint,  which  is  also  an 
assumption,  is  that  when  wq  *  w^/2,  A  must  equal  B.  From 
the  first  constraint,  this  implies  that  for  Wq  « 

A  =  B  »  0.5  (15) 

The  magnitudes  of  X  from  Figure  9  are 

Y  -  X  (Wo)  »  M  *  slnii(Wn/Ws)  (16) 

'(Wq/Ws) 

and 


Z  -  X  (Ws-Wq)  »  M  *  sinic(ws  -Wo)/Ws  (17) 

K  (  Ws  -Wq  )  /  Ws 
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Assuming  that  A  and  B  are  similar  sine  functions,  where 


sinc(u) 


(8in(u)  ]/u,  define 


To  use  Figure  17,  a  sample  rate  and  Input  freq¬ 
uency  (for  a  sine  or  cosine  input)  is  selected  and  the  ratio 


(i>q/(i>s  is  found.  Then,  from  Figure  17,  B  ,  the  magnitude 
of  the  impostor  frequency  is  found.  Rearranging  Eq.  (14) 

A  =  1-B  (20) 

the  new  magnitude  of  the  input  frequency.  A,  attenuated  by 
the  ZOH  is  calculated. 


Model  Examples 


Example 

One .  Let 

Ws  = 

2»k*64 

(21 

II 

2*11*31 

(22 

T  * 

2*ii/Ws  »  0.015625 

for  an  input 

r(t)  = 

M*cos (Uo*t ) 

with  M=l,  the 

sampled 

output,  from  MATRIX  , 

shown  in  Figure 

18(a),  is 

c(kT)  = 

cos  ( 2*ii*31*)c*T ) 

(  23 

The  model,  from  equation  (13),  is 

c^(t)  =  A*cos (Wot ) +B*cos ( Ws -Wq ) t 


where 


T  (SECONDS) 


3? 


From  Figure  17,  B=.48,  so  that  A=l-.48=.52.  The  model  is 
now  defined  as 

c^(t)  =  .52*cos(2»ii*31»t)  +  .48*co8(2»ii*33*t)  (25 

Figure  18(b)  shows  the  model  output. 

Example  Two .  Let 


Us  =  2*k»64  (26 

Uq  =  2*k*21  (27 

T  =  2»ii/Ws  =  0.015625 

The  sampled  output,  shown  in  Figures  19(a)  and  20(a),  is 

c(kT)  =  cos( 2*n*21*k»T)  (29 

The  model  from  equation  (13)  and  Figure  17 
for 

«s  =  2*n*64  -  2*n*21  =  43 

A  =  .77 

B  =  .23 

Is  now  defined  as 

Cm(t)  =  .77*cos(2*ii*21*t)  +  .23*cos(2*x*43*t)  (30) 


Figures  19(b)  and  20(b)  show  the  model  output. 


T  (SECONDS) 


Figure  19(a).  Output  c(kt)  =  cos( 2*ii»21»k*T)  Subharmonic 
Allas  from  Sampled  Signal  for  £o=21  Hz 
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Figure  19(b).  c(t)=  0 . 77*cos  (  2*ii*21*t ) +0 . 23*cos  (  2*ii*43*t ) 

Subharmonic  Aliasing  Model  for  fQ=21  Hz 
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Figure  20(a).  Output  c(kt)  *  cos(2*ii»21»k*T)  Over  several 
Subharmonic  Alias  Cycles  for  fo®21  Hz 


Three .  Let 


Wg  =  2*11*1  (31) 

«o  =  2*k*0.21  (32) 

M  =  5 

for  an  input 

r(t)  =  5*co8(2*x*0.21*t) 

The  sampled  output,  shown  in  Figure  21(a),  is  given  as 

c()tT)  =  5*cos(2*n*0.21*k*T)  (33) 


The  model,  from  equation  13,  is  given  by 


Cj^(t)  =  5*[A*cos(«Qt)+B*cos(«s-«Q  )t ) 

where 

Ws-Wq  «  2*K*(1-0.21) 

=  2*x*0.79 
Wq/Us  ®  0.21 


From  Figure  17,  B=.l,  so  that  A=l-.1=.9.  The  model 
is  now  defined  as 

c^(t)  =  5*l0.9*cos(2*n*0.21*t) 

+  .l*cos(2»x»0.79*t)  J  (34) 


Figure  21(b)  shows  the  model  output. 

Having  validated  the  model,  it  can  now  be  used  to 
determine  the  frequency  and  magnitude  of  the  subharmonic 
aliased  signal.  Two  phenomena  occur  when  two  sinusoids  are 
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added  together,  they  are  beating  and  apparent  low  frequency 
surge . 

Beating  occurs  when  the  ratio  of  the  two  sinusoidal 
frequencies  Is  nearly  unity  (13:21).  Example  1,  used 
previously,  shows  this  phenomena  well.  The  model 
frequencies  are  given  as 

fo  =  31  Hz 

fg-fQ  =  33  Hz  (36) 

The  beat  frequency  f^  is  the  absolute  magnitude  of  the 
difference  between  these  two  frequencies. 

f.  =  33-31  =  2  Hz  (37) 

D 

The  amplitude  of  the  signal  varies  between  the  sum  and 
the  difference  of  the  component  amplitudes  where 

M*A  =  1»0.52  =  0.52  (38) 

M*B  =  1*0.48  =  0.48  (39) 

so  that  the  amplitude  varies  between 

0.52  +  0.48  =  1  (40) 

and 

0.52  -  0.48  =  0.04  (41) 

So,  the  resultant  wave  has  the  same  apparent  frequency 
as  the  major  component  (ie.  that  with  the  greater  amplitude) 
(13:21)  of  31  Hz,  its  amplitude  varies  between  1.0  and  .04, 


and  It  has  a  beat  frequency  of  2  Hz.  This  Is  shown  In 
Figure  22. 


Figure  22.  c(t)=  0 . 52*cos  (  2*TC*31»t ) +0 . 48*cos  (  2*ii*33*t ) 

Subharmonlc  Allas  with  a  Beat  Frequency  for  £(,  =  31  Hz 


An  apparent  low-frequency  surge  occurs  when  the  ratio 
of  the  two  sinusoidal  frequencies  Is  nearly  an  Interger  (n) 
(when  the  larger  frequency  Is  divided  by  the  smaller 
frequency)  (13:21).  Example  2,  used  previously,  shows  this 
phenomena  well.  The  model  frequencies  are 


21 

Hz 

(42 

43 

Hz 

(43 

The  apparent  low-frequency  surge  f-i  Is  the  absolute 


magnitude  of  the  difference  between  n»fQ  and  fg-f©  (13:39), 
where  n.  In  this  case  Is  equal  to  two.  Therefore,  the 
apparent  low-frequency  surge  Is 


fg,  =  43-(2*21)  (44) 

=  1  Hz 

Again,  the  amplitude  of  the  signal  varies  between  ♦’he 
sum  and  the  difference  of  the  component  amplitudes  where 

M»A  =  1*0.77  =  0.77  (45) 

M*B  =  1*0.23  =  0.23  (46) 

so  that  the  amplitude  varies  between 

0.77  +  0.23  =  1  (47) 

and 

0.77  -  0.23  =  0.54  (48) 

So,  the  resultant  wave  has  the  same  apparent  frequency 
as  the  major  component  (ie.  that  with  the  greater  amplitude) 
of  21  Hz,  its  amplitude  varies  between  1.0  and  .54,  and  it 
has  an  apparent  low-frequency  surge  of  1  Hz.  This  is  shown 
in  Figure  23. 
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Figure  23.  Subharmonlc  Aliasing  with  an  Apparent  Low 

Frequency  Surge 


Utility  of  the  Model 

Now  that  the  model  is  well  defined,  how  can  it  be  used? 
One  possible  application  Is  to  use  the  subharmonlc  alias 
model  as  an  Input  (known  to  be  within  the  subharmonlc 
range,  <  Wq  ^  so  that  the  output  characteristics 

of  the  system  can  be  predicted.  A  simple  system  with  a 
first  order  filter  Is  used  to  demonstrate  this. 

First  Order  Filter  Example .  The  purpose  of  this 
example  Is  to  show  that  the  subharmonlc  aliasing  model  can 


be  used  much  like  the  pseudo  continuous  time  model  (8:246) 
to  predict  the  output  characteristics  of  a  digital  system. 

First/  the  system  In  Figure  24(a)  Is  simulated  on  a 
hybrid  computer  to  determine  the  base-line  system  response, 
that  Is,  the  actual  response  of  the  anolog/dlgltal  control 
system.  The  system  parameters  are 

r(t)  -  cos(2»ii*31*t)  volts 
T  =  0.015625  seconds 
Mg  =  2*ii/T  =  2*ii»64  (radians/second) 

where  Wq  =  2»ii*31  of  r(t)  is  within  the  subharmonic  range, 
2*11*6.4  <  Uq  <  2*11*32,  for  this  example.  The  steady  state 
output  response  for  this  system  is  found  by  determining  the 
frequency  response  of  the  first  order  filter.  For  the  input 
frequency  Wq  »  2*ii*31  the  output  Is  attenuated  by  »  -7  dB 
with  a  phase  shift  ^  =  -63°  .  Therefore,  the  output  equation 
for  this  system  Is 

C(kT)  =  .45*cos(2*it*31*t  -  63°)  (49) 

This  output,  shown  in  Figure  25  and  labeled  C(kT),  is  the 
hybrid  computer  response  to  the  system  modeled  in 
Figure  24(a)  and  correlates  well  with  output  equation  (49). 
Unfortunately,  equation  (49)  says  very  little  about  the 
subharmonic  characteristics  of  the  output  In  Figure  25  which 
definitely  shows  a  beat  frequency  of  2  Hz  with  an  amplitude 
varying  from  approximately  0.02  to  o.45  volts. 
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Figure  25.  Outputs  from  Hybrid  Simulation  C(kT),  Pseudo 
Continuous  Time  Simulation  CpQj,  and  Model  Simulation  C|^ 


50 


approximated  by  using  1/T  and  the  first  order  Fade' 
approximation  (8:249)  respectively  to  form 

Ga  =  2/T  (50) 

s  +2/T 

Substituting  the  sample  period  T=0. 015625  seconds  into  (50) 

G.  =  128  (51) 

s  +  128 

The  steady  state  system  output  response  Cp^j  for  the  pseudo 
continuous  time  second  order  system  model  is 

*  0.249*cos(  2*it*31*t  -  120°)  (52) 

This  output  response  is  shown  in  Figure  25  and  is  labeled 
CpcT*  This  is  not  a  good  model  of  this  system  for  two 
reasons.  First,  it  does  not  reflect  the  existence  of  a 
subharmonic  alias,  which  is  linown  to  exist,  as  shown  in 
output  C(kt).  Second,  the  model  itself  lies  outside  the 
"good"  Tustln  region ( 8 : 249 ,  which  is  described  as  the 
s-domain  region  which  adequately  models  a  discrete  system 
in  the  continuous  domain  so  that  continuous  control  law 
can  be  used  to  determine  the  system  response. 

Third,  model  the  input  using  the  subharmonic  aliasing 
aliasing  model 


r|^(t)  =  A*cos(Wo^)  *’  B»co8  («s -Wq  )  t  (53) 


Wfir  I 


.  •  ^  • 


.  •  «  •  (5  4 


Wq  Ws  =  l  >  4 


SO  tro»  Figure  1' 


B  =  0.48 


and  A  =  1 -B  =  0.52 


The  Input  r(t)  la  now  defined  as 


r(ti  =  0 . 52*cos(  2*ii*31*t ) 


+  0.48*cos(2*K*33*t) 


For  linear  systems  (which  this  is)  the  principle  of 
superposition  can  be  used  to  calculate  the  output.  The  two 
components  of  the  modelled  input  can  be  considered  one  at  a 
time.  The  output  signal  is  then  equal  to  the  sum  of  the 
contributions  produced  by  each  input  (15:162).  Let  Cj^-j  be 
the  output  response  to  0 . 52*cos  (  2*n*31*t )  and  C|^2 
output  response  to  0 . 48*cos ( 2*n*33*t ) .  Therefore,  the 
steady  state  responses  are 


C,^.,  =  .235»cos(2»ii*31»t  -  63  ) 
C,^2  =  .206*cos(2*n*33»t  -  64°) 


Combining  (55)  and  (56)  the  output  CM  becomes 


C,^(t)  *  0.235*cos(2»it»31*t  -  63°) 
+  0.206»cos(2»ii*33*t  -  64°) 


Ml  i)a 


‘  w_ 


This  output,  shown  In  Figure  25  and  labeled  as  shows 

a  definite  subharmonic  alias  whose  characteristics  can  now 
be  calculated.  From  (57)  the  beat  frequency  can  be  calcu¬ 
lated  as  (33-31)  Hz  =  2  Hz.  The  amplitude  of  the  beat 
frequency  varies  between  0.441  and  0.029. 

This  example  demonstrates  that  the  subharmonic  aliasing 
input  model  produces  the  salient  characteristics  of  the 
digital  (or  hybrid)  subharmonic  aliased  output.  The  model 
output  can  then  be  used  to  predict  the  subharmonic  aliasing 
character ist ics . 
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Interchannel  Difference 


The  purpose  of  this  chapter  Is  to  provides  a  means  to 
predict  the  Interchannel  difference  for  Input  frequencies 
within  the  subharmonic  range.  In  Chapter  II,  the  maximum 
sample  period  Is  calculated  with  Equation  (1)  so  that  the 
interchannel  difference  would  be  less  than  five  percent 

where  T  is  the  maximum  sample  period,  Is  the  maximum 

interchannel  difference,  and  s  is  the  maximum  signal  rate  of 
change.  As  described  in  the  Data  Aquisitlon  section  of 
Chapter  II,  the  sensor  sample  rate  T  for  the  AFTI/F-16  Is 
0.00390625  seconds.  Since  the  sensor  sample  rate  (1/T)  Is 
fixed,  the  maximum  Interchannel  difference  for  an  Input  can 
easily  be  calculated  by  rearranging  Equation  (1)  to 

emax  *  T  *  s 

Given  a  sinusoidal  input 

r(t)  =  M*cos(«(j*t) 

the  maximum  rate  of  change  s  is  found  by  differentiating 

r(t)  with  respect  to  time 

s  =  maxidr (t)/dt}  (60) 

Wot 

s  =  max  {M*WQ*sin(WQ*t) } 

Wot 

=  +  M  *  Wq 

for  Uot=  ±nn/2  n=l,3,5,... 


(58) 

(59) 


P 


-.-"V 


- 

S;v, 


***** 


■/ 
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Substituting  Equation  (59)  into  Equation  (58)  for  n^l  gives 


e  s  T  ^  M 
®max  *  o 


Let  e  be  defined  as  the  maximum  Interchannel  difference  when 


M=l,  therefore. 


e  =  T  •  «r 


Substituting  (62)  Into  (61) 


emax  =  «  *  e 


A  graph  of  Equation  (62)  shows  the  relationship  between  the 
frequency  of  an  Input  signal  (Wq)  and  the  maximum  error  e^gj^ 
(normalized)  caused  by  time  skev  between  computers  for  the 
sensor  sample  period  of  3.90625  msec  (Figure  26).  Prom  this 
figure.  It  Is  obvious  that  the  high  frequencies  (le.  freq¬ 
uencies  In  the  subharmonic  range)  cause  large  Interchannel 
differences,  which  may  exceed  the  Input  Selector  Monitor 
trip  levels.  If  this  frequency  (Input  signal)  persists 
sufficiently  long,  a  sensor  may  be  voted  as  failed  and 
removed  from  the  failure  management  scheme. 


Example 


The  Inter-branch  trip  level  Is  15%  of  full  scale  for 
most  AFTI/F-16  inputs.  This  example  uses  a  hypothetical 
Input  of  -4.5  to  4.5  volts  full  scale  and  a  trip  level  (TL) 


'“O  *  ' 


NORMALIZED  ERROR 


Figure  26,  Maximum  Interchannel  Difference 
for  Maximum  Time  Skew( . 00390625  seconds) 


of  +  15  percent,  so  that  the  trip  level  for  this  example  is 

TL  =  +  9  »  0,15  (64) 

=  +  1,35 


Given  an  input 


r(t)  =  3*cos(2*ii»21,186*t)  (volts) 
The  maximum  rate  of  change  for  this  input  is 

s  =  3*2*tc21.186  (volts) 


The  sensor  sample  period  T  Is  0.00390625  seconds,  therefore 
the  maximum  Interchannel  difference  from  Equation  (56)  is 


) 

«max  '  -  0.00390625»3»2»«*21.1816 
«  +  1.56  volts 


which  is  greater  than  the  trip  level  of  +.  1.35  volts.  This 
implies  that  for  this  input  magnitude  (M  «  3)  and  frequency 
(Uq  =  21.186  Hz)  if  two  FLCC's  are  time  s)iewed  by  the 
maximum  amount  (0.00390625  seconds)  a  transient  error  occurs 
at  least  twice  per  cycle.  This  example  is  simulated  on  a 
digital  computer  MATRIX  and  a  hybrid  computer  (SIHSTAR 
Programs,  SI. SATL, Appendix  B)  for  maximum  time  sicew  with 
results  shown  in  Figure  27  and  Table  III.  Figure  27(a)  is 
the  time  response  of  the  input  r(t)  =  3»cos ( 2*x*21 . 186 ) . 
Figure  27(b)  is  the  time  response  of  the  Interchannel 
difference  for  maximum  time  skew.  This  figure  is  generated 
by  sampling  the  input  r(t)  every  0.00390625  seconds  and 
subtracting  the  roost  recent  past  sampled  value  by  the 
present  sampled  value.  Figure  27(b)  clearly  demonstrates 
that  the  trip  level  (drawn  in  at  +  1.35)  will  be  exceeded  at 
least  eight  times  within  the  first  0.1  seconds  for  this 
input.  The  data  items  marked  by  three  asterisks  in  Table 
III,  show  that  there  are  often  more  than  two  trips  per  input 
cycle.  This  chapter  demonstrates  the  sensitivity  of  the 
interchannel  difference  with  respect  to  frequencies  in  the 
subharmonic  range  and  the  ease  of  calculating  the  maximum 


27(a).  Input  Frequency  (21.186  Hz)  Sampled 
625  sec.  (the  Sensor  Sample  Period  for  the 


Table  III.  Interchannel  Difference  for  Maximum  Time  Skew 
(Input  r(t)=3*cos(2»ii»21.186)  sampled  every . 00390625  sec.) 


LIME 

TIME 

SAMPLED 

INPUT 

INTERCHANNEL 

DIFFERENCE 

0 

o 

• 

o 

0.00000 

0.0 

1 

0.003906 

2.20647 

-2.20647 

*** 

2 

0.007813 

2.9983 

-0.791831 

3 

0.011719 

2.13425 

0.864045 

4 

0.015625 

0.094183 

2.04007 

*** 

5 

0.019513 

-1.99803 

2.09221 

*** 

6 

0.023438 

-2.98926 

0.991235 

7 

0.027344 

-2.32777 

-0.661489 

8 

0.03125 

-0.3062 

-2.02157 

*** 

9 

0.035156 

1.78601 

-2.09221 

*** 

10 

0.039063 

2.95152 

-1.16551 

11 

0.042969 

2.50258 

0.448946 

12 

0.046875 

0.672736 

1.82984 

*«* 

13 

0.050781 

-1.52711 

2.19985 

14 

0.054688 

-2.88436 

1.35724 

*** 

15 


0.058594 


-2.65068 


-0.233671 


The  purpose  of  this  chapter  Is  to  demonstrate  the  rate 
limiting  effect  that  the  AFTI/F-16  engineers  described  as 
"the  rate  limiter  following  the  low  frequency  alias 
waveforms,  instead  of  the  input  wavef orm" . ( 1 : 10 )  Also,  for 
the  benefit  of  AFTI/F-16  engineers,  this  chapter  provides 
the  output  for  comparison  of  two  different  rate  limiter 
configurations:  the  inputs  averaged  then  rate  limited,  the 
Inputs  rate  limited  then  averaged.  Hybrid  computer  simula¬ 
tions  are  used  to  provide  the  outputs. 

Chapter  IV  shows  how  the  analog  to  digital  and  digital 
to  analog  processes  induce  impostor  frequencies  which 
distort  the  original  input.  At  predetermined  frequencies 
(based  on  the  sample  rate)  well  defined  subharroonic  aliases 
are  formed,  with  beat  frequencies  or  apparent  low  frequency 
surges.  Chapter  V  verifies  that  inputs  in  the  subharmonic 
range  ^  “5/2  have  a  high  rate  of  change,  which 

effects  the  magnitude  of  the  interchannel  difference. 

The  purpose  of  the  digital  flight  control  system  soft¬ 
ware  rate  limiter  is  to  limit  this  rate  of  change  to  no 
more  than  five  percent  full  scale  per  computational  frame. 
The  combination  of  three  elements: 

1.  the  apparent  low  frequency  surge, 

2.  the  input  high  rate  of  change  and. 


3.  the  software  rate  limiters  nonlinear  characteristics 

produce  the  effect  of  the  rate  limiter  following  the 
"low  frequency  component". 

Since  the  DFCS  control  laws  and  aircraft  dynamics  are 
not  simulated,  a  generic  input  to  the  rate  limiter  is  used 
to  simulate  the  effect  of  the  rate  limiter.  To  determine 
a  reasonable  rate  limit,  the  example  from  Chapter  V  is  used 
again.  Full  scale  is  defined  as  -4.5  to  4.5  volts.  To 
limit  this  to  five  percent  per  computation  frame  the  rate 
limit  threshold 

Tpi^  =  9  volts  *0.05 
=  0.45  volts 

The  pseudo-code  for  rate  limiter  is  provided  below 
and  is  also  shown  in  block  diagram  form  in  Figure  28. 

if  old_output  -  new_lnput  >  rate_llmlt_threshold  then 
new_output  *  old_output  -  rate_limit_threshold 
else  if  old_output  -  new_lnput  <  -  rate_llrolt_threshold  then 
new_output  »  old_output  rate_l  lrolt_threshold 
else  new_output  =  new_lnput 
end  if 

old_output  »  new_output 

where  old_output,  (Shown  in  Figure  28  as  Ol(k-l)*T)l)  is  the 
rate  limiters  previously  computed  output,  new_output 
(0{k*T)>  is  the  the  output  of  the  rate  limiter,  new_input 


Figure  28.  Software  Rate  Limiter  Block  Diagram 

{I(k*T)i  is  the  input  to  the  rate  limiter,  4(k»T)  is  the 
difference  between  old_output  and  new_input,  and  is  the 

rate_llmlt_threshold.  This  pseudo-code  is  similar  to  the 
FORTRAN  code  used  in  the  hybrid  simulations. 

Three  software  rate  limiting  simulations  are  performed 
on  a  hybrid  computer  (SIMSTAR  Programs,  Sl.RLSIM,  Appendix 
B)  to  demonstrate  the  rate  limiting  effect  on  the  subharmonic 
alias.  The  first  simulation  (Block  Diagram,  Figure  29(a)) 
rate  limits  a  single  input  simply  to  show  the  rate  limiter 
following  the  "low  frequency"  of  the  subharmonic  aliased 
Input).  The  second  simulation  (Figure  29(b))  provide  two 
inputs,  time  skewed  by  one  sample  period  (0.015625  second, 
which  is  approximately  one  FLCC  frame)  which  are  averaged, 
then  input  to  the  rate  limiter.  The  third  simulation 
(Figure  29(c))  rate  limited  the  two  time  skewed  signals  and 
then  averaged  them. 


To  Illustrate  the  effects  of  the  rate  limiter,  a  low 
frequency  is  sampled  at  a  low  sample  rate  to  induce 
subharmonic  aliases  and,  at  the  same  time,  minimize  the 
strip  chart  recorder  pin  dynamics.  Figure  30  sho%<s  the 
input,  D/A  output,  and  the  rate  limiter  response.  The  input 

r(t)  ■  co8(2*ii*fo»t) 

where  f^  =  1  Hz  sampled  at  3.0476  samples/sec.  From  Chapter 
4,  this  input  and  sample  rate  produce  a  well  defined  subhar¬ 
monic  alias  with  an  apparent  low  frequency  surge  of 

fo  *  (3*1)  -  3.0476 
=.0476  Hz 

The  rate  limit  threshold  is  set  at  0.45  volts  per 
computation  frame,  which,  for  a  step  input  with  magnitude  2 
to  the  rate  limiter  requires  5  samples  for  the  output  of  the 
rate  limiter  to  reach  2.  This  is  calculate  by  dividing  the 
magnitude  of  the  step  input  by  the  rate  limit  threshold  and 
rounding  up  to  the  nearest  integer,  for  example,  2/0.45  » 
4.44,  rounding  up  =  5.  With  the  aid  of  the  rate  limiter 
pseudo-code  and  points  labeled  on  Figure  34  the  effect  of 
rate  limit  is  described.  In  general,  for  a  rate  limit 
threshold  of  Tp^^  =  .45,  comparing  old_outputs  (A'  -  E'  in 
Figure  30)  with  new  inputs  (B  -  E  in  Figure  30)  the 


Input  (iHz)  U 


f- 


Subharmonlc  0 
Alias 


Rate  Limiter 
Output 


\re' 


U'  C' 


4.0  6.0 

t  ( seconds ) 


Figure  30.  Rate  Limiter  Output  to  1  Hz  Input 
Sampled  at  3.0476  samples/second 
(A,  B,  C,  D,  and  B  are  values  of  the  subharmonic  alias) 
(A',B',  C,  D',  are  E'  are  rate  limiter  output  values) 
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new_output8  can  be  calculated.  For  this  example  A*  -  B  is 
less  than  the  rate_l lrolt_threshold 

A'  -  B  <  -T 

so  that  the  new_output  is 

B’  -  A'  +  T 

Simi lar ly^ 

B'  -  C  >  T 

80  that  the  new_output  is 

C*  «  B*  -T 

and 

C  -  D  <  T 

so  that  the  new_output  is  equal  to  the  new_lnput 

D'  «  D 

These  three  rate  limits  are  examples  of  the  three 


Figure  31  is  this  example  run  on  a  hybrid  computer 
(SIMSTAR)  with  a  slower  strip  chart  recorder  speed  to  show 
two  full  cycles  of  both  the  subharmonic  alias  and  the  rate 
limiter  output.  This  is  an  good  example  of  how  the  non¬ 
linear  characteristics  of  the  software  rate  limiter,  and 
the  characteristics  of  the  subharmonic  alias  combine  so  that 
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the  rate  limiter  follows  the  "low  frequency". 
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Compar Ison  o£  Techniques 

This  section  compares  the  outputs  of  the  two  different 
possible  rate  limiting  techniques:  averaging  the  Inputs  and 
then  rate  limiting  (C  ,  Figure  29(b)),  limiting  the  inputs 
and  then  averaging  (C  ,  Figure  29(c)).  Inputs  are  of  the 
form  r(t)  »  cos((iiQ*t),  sampled  at  64  samples/second  with 
with  a  rate  limit  threshold  Tp|^  =  0.45. 

For  Wq  *2*k*31  (Figure  32),  both  techniques  effectively 
rate  limit  the  Input  signal.  For  Wq  *2***2!  (Flgure33), 
rate  limiting  the  Inputs  prior  to  averaging  better  limits 
the  high  frequency,  but  the  magnitude  of  the  low  frequency 
is  greater  than  that  of  the  configuration  which  averages  the 
Inputs  prior  to  rate  limiting.  For  Wq  =2*ii*17  (Figure  34), 
rate  limiting  prior  to  averaging  decreases  the  magnitude  of 
both  the  high  and  low  frequencies  better  than  averaging  and 
then  rate  limiting.  For  Wq  =2*ii*13  Hz,  another  well 
defined  subharmonic  alias  (Figure  35),  rate  limiting  prior 
to  averaging,  once  aglan,  better  limits  the  Input  signals. 
Finally,  at  .5  Hz  (Figure  36),  (which  Is  near  the  short 
period  mode  of  the  AFTI/F-16)  there  is  no  limiting  for  this 
Input  -  which  Is  expected  considering  the  low  rate  of  change 
of  the  Input. 

In  general,  rate  limiting  the  Inputs  prior  to  averaging 
limits  the  Input  signal  better  than  averaging  the  Input  and 
then  rate  limiting. 
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Figure  32.  Rate  Limiter  Output  to  31  Hz  Input 
Sampled  at  64  saroples/sec 
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Figure  35.  Rate  Limiter  Output  to  13  Hz  Input 
Sampled  at  64  samples/second 


VI I .  Conclusions  and  Reconunendat Ions 

Conclusions 

There  are  five  main  points  covered  In  this  thesis:  first, 
determine  the  cause  of  subharmonic  aliasing,  as  described  by 
the  AFTI/F-16  engineers  as  "the  creation  of  uncorrelated 
low  frequencies  whenever  a  subharmonic  of  the  sample 
frequency  Is  Input  Into  the  system";  second,  model  the 
subharmonic  aliasing  effect,  so  that,  by  knowing  only  the 
Input  frequency  and  the  system  sample  rate  the  output 
characteristics  can  be  calculated;  third,  demonstrate  the 
utility  of  the  subharmonic  aliasing  model  when  used  as  an 
Input;  fourth,  mathematically  show  how  Inputs  In  the 
subharmonic  range  (Wg/10  <  «q  <  may  Induce  large 

Interchannel  differences;  and  fifth,  demonstrate  the  effect 
of  the  rate  limiter  on  the  subharmonic  alias  and  compare  two 
rate  limiter  configurations. 

The  mathematical  model  determined  that  subharmonic 
aliasing  Is  the  result  of  Impostor  frequencies  (Wg-  Wq, 
much  like  aliasing)  being  Introduced  Into  the  output  signal 
by  the  sampling  process.  The  magnitude  of  the  Impostor  Is 
dependant  upon  the  characteristics  of  Input  filters  and  the 


ZOH.  Well  defined  subharmonic  aliases  occur  due  to;  the 
impostor  frequency  and  one  of  two  phenomena,  either  the 
apparent  low  frequency  surge  or  the  beat  frequency.  The 
apparent  low  frequency  surge  occurs  when  the  the  input 
frequency  is  nearly  an  Integer  multiple  of  the  sampling 
frequency.  The  beat  frequency  occurs  when  the  input 
frequency  is  nearly  one  half  of  the  sampling  frequency. 

Simulations  from  Chapter  IV  demonstrated  the 
utility  of  the  model  as  an  input  used  to  predict  the 
characteristics  of  the  subharmonic  aliased  signal. 

Mathematical  calculations  and  simulations  in  Chapter  V 
verify  that  for  even  small  magnitude  inputs  large  inter- 
channel  differences  may  occur  for  frequencies  in  the  subhar¬ 
monic  range. 

Simulations  discussed  in  Chapter  VI  show  that  it  is  the 
combination  of  well  defined  subharmonic  aliases,  the 
relative  high  frequencies  of  the  signal,  and  the  nonlinear 
characteristics  of  the  rate  limiter  which  form  the  low 
frequency  output  from  the  rate  limiter.  Simulations  also 
that  rate  limiting  input  signals  prior  to  averaging 
them  provid  a  better  overall  rate  limiting  effect. 


Recommendations 


Since  subharmonic  aliasing  occurs  for  high  frequencies 
(relative  to  the  sample  rate),  sampling  at  a  higher  rate 
reduces  the  undesirable  subharmonic  aliasing  effect. 
Additionally,  since  the  phenomenon  Is  actually  the 
combination  of  two  high  frequencies,  using  higher  order 
analog  filters  and/or  filters  with  lower  cutoff  frequencies 
will  attenuate  the  Input  signals  in  this  range. 

Unfortunately,  the  two  simple  solutions  recommended 
above  are  not  easily  implemented,  since  It  Is  very  difficult 
to  change  the  sample  and  frame  rates  once  a  computer 
architecture  Is  fixed,  a^id  since  higher  order  (and/or  lower 
cutoff  frequency)  filters  tend  to  reduce  the  gain  and  phase 
margins  of  the  control  system. 

Further  research  should  involve  simulating  combinations 
of  the  sensor  sample  rates.  Input  and  Output  selector 
monitor  thresholds  limits.  Rate  Limiter  configurations, 
anti-aliasing  filter  cutoff  frequencies,  and  computation 
frame  rates  to  determine  optimal  thresholds  for  the  system. 
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Appendix  A;  Examples  of  Subharmonic  Aliasing 

Figures  A1-A7  are  examples  of  Subharmonic  Aliasing  with 
either  a  beat  freqeuncy  or  an  apparent  low  frequency  surge. 
These  plots  are  generated  with  an  Input  r  ( t )  =cos  ( 2»Ti*fQ*t ) , 
where  f^  Is  the  Input  frequency  In  Hz,  to  an  analog  to 
digital  converter  and  then  to  a  digital  to  analog  converter 


1  sec 


Figure  A1 .  Subharmonic  Aliasing  Effect  (input  29  Hz, 

sampled  at  64  samples/sec) 


Figure  A2 .  Subharmonic  Aliasing  Effect  (Input  26  Hz, 

sampled  at  64  samples/sec) 


Figure  A3.  Subharmonic  Aliasing  Effect  (Input  20  Hz, 

sampled  at  64  samples/sec) 
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Figure  A5. 
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Subharmonlc  Aliasing  Effect  (input  15  Hz, 
sampled  at  64  samples/sec) 
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Figure  A6 ,  Subharmonic  Aliasing  Effect  (input  11  Hz 

sampled  at  64  samples/sec) 
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Figure  A7 .  Subharmonic  Aliasing  Effect  (Input  7  Hz, 

sampled  at  64  samples/sec) 


Figures  A8-A13  are  examples  of  a  frequency  which  is  a 
multiple  (1/3  in  thl-j  case)  of  the  sample  frequency.  These 
figures  demonstrate  how  the  phase  If  the  input  determines 
the  magnitude  of  the  output. 


Figure  A8 .  Subharmonic  Aliasing  Effect  (input  21.333  Hz, 
0°  phase  shift,  sampled  at  64  samples/sec) 
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Figure  AlO.  Subharmonic  Aliasing  Effect  (input  21.333  Hz 
20°  phase  shift,  sampled  at  64  samples/sec) 
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Figure  All.  Subharmonic  Aliasing  Effect  (input  21.333  Hz 
30°  phased  sh'ft,  sampled  at  64  samples/sec) 
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Figure  A12.  Subharmonic  Aliasing  Effect  (input  21.333  Hz 
40°  phased  shift,  sampled  at  64  samples/sec) 
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Appendix  B;  SIMSTAR  PROGRAMS 

The  following  programs  were  used  to  demonstrate 
subharmonic  aliasing,  and  Its  effects  using  the  SIMSTSAR 
hybrid  computer.  Sl.SASIM  (Page  88)  simulates  subharmonic 
aliasing  with  a  simple  A/D  -  D/A  circuit  sampling  a  sinu¬ 
soidal  input.  Sl.SAPIL  (Page  91)  simulates  the  effect  of 
an  input  filter  (FILTER  in  radians)  on  the  subharmonic 
alias.  SI. TESTS  (Page  93)  simulates  the  effect  of  an 
output  filter  (FILTER  in  radians)  on  the  subharmonic  alias. 
Sl.SATL  (Page  95)  simulates  the  maximum  interchannel 
difference  due  to  time  slcew  for  user  specified  threshold 
levels  (TL).  Sl.RLSIM  (Page  97)  simulates  the  effect  of  a 
software  rate  limiter  on  the  subharmonic  alias  for  three 
specific  cases:  a  single  input  rate  limited,  two  inputs 
(the  past  input  and  the  present  input)  averaged  together 
and  then  rate  limited,  and  the  inputs  rate  limited  and  then 
averaged  together. 
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*PSP=1,0,ERR=ALL 

♦TITLE 

THESIS  SUBHARMONIC  ALIASING  SIMULATION 
♦INPUT 

'This  Is  a  self  contained  program  and  requires  no  ' 
'external  inputs.  All  outputs  go  to  the  strip' 
'chart  recorder' 

I  I 
I  I 

'If  you  are  typing  this  program  in,  do  not  include 
'line  numbers' 

'THIS  PROGRAMS  INPUT  IS  OF  R=SIN ( 2^PI ♦FREQ^T  +  PHI)' 

' WHERE ' 

'  FREQ  =  THE  INPUT  FREQUENCY  IN  Hz  AND  IS  SET' 

'  BY  THE  OPERATOR' 

'  PHI  =  THE  PHASE  SHIFT  IN  DEGREES  AND  IS  SET' 

'  BY  THE  OPERATOR' 

'  THE  SAMPLE  RATE  IS  64  SAMPLES/SECOND  AMD  MAY  BE' 

' CHANGED  BY  THE  OPERATOR . ' 

'  TO  DO  THIS' 

•SET  CINT=sample  period  required' 

'SET  PERIOD=sample  period  required' 

PROGRAM 

' INTERRUPT  DECLARATIONS ' 

' INTDEF(0,1,1) ' 

'INTDEF(1,1,0) ' 

INITIAL 

'eBETA(BETA) ' 

MAXT  *  PERIOD/BETA 
LOGPER  =  Cl NT  ♦  BETA 
FRAD=FREQ^2^PI 
PRAD=PHI^PI/180 
VARIABLE  TIME  =  0 

CONSTANT  BETA=1  ,  RUNTIM  =  3.0  ,  PERIOD=  .015625 
CONSTANT  CINT=. 015625,  FREQ=32,  PHI =0, PI =3 . 14159 
'ePARAMETER  BETA, RUNTIM, FREQ, PHI ' 

'§MAXVAL  FREQ=  64,PHI=360' 

'iMINVAL  FREQ=.01,PHI=0' 

'eMAXVAL  BETA  =100,  RUNTIM  =50,  TIME=50' 

'gMINVAL  BETA=.l  ,  RUNTIM  =0  ,  TIME=0  ' 

NSTEPS  NSTEP=1 
D1  =  0 

END  $'OF  INITIAL' 

DYNAMIC 

'INTERRUPT  RATE  ERROR  DECLARATIONS' 

LOGICAL  ENDERl,  RATERl,  ERRORl 
ENDERl  =  .FALSE. 

ERRORl  =  RATERl 
DERIVATIVE 

CONSTANT  D1MAX=1 
CONSTANT  A1MAX=1 
CONSTANT  CMAX=1 
'§SCALE  D1=D1MAX' 


'eSCALE  Al^AlMAX* 

•eSCALE  C=CMAX' 

D1=A1 

'^PARALLEL' 

R’’HARM(0,FRAD,PRAD) 

A1*R 

C=D1 

TERMT(TIME  .GT.  RUNTIM) 

• ©RECORD (REC01,,,,,,,,R,C,A1,D1,,,, ) ' 
•DEFINE  INTERRUPT  CONTROL' 

GPIO  =  CLOCK (PERIOD) 

GPU  =  CLOCK  (LOGPER) 

'©INTRRT  1  =  GPIO' 

•eiNTRRT  2  =  GPU' 

RATBRl =RATERR (GPIO, ENDERl ) 

'©END  PARALLEL' 

END  $'OF  DERIVATIVE' 

END  $'OF  DYNAMIC 
TERMINAL 

END  $'OF  TERMINAL' 

END  $'OF  PROGRAM' 

*TRANSLATE 

DCA(1)=D1 

PADC(1)=A1 

♦OUTPUT 

*END 

SUBROUTINE  PREPl 

+ 

INCLUDE  El.SASIM 
A1*QRPADC(0)*S;A1 
RETURN 
END 
C 

SUBROUTINE  POSTl 

+ 

INCLUDE  El.SASIM 
COMMON  /QQDCP/DCASF (0:0) 

LOGICAL  DELAY 

CALL  QWDCAR(0,D1*DCASF(0) ) 

IF  (L:RATER1)  CALL  ZZRTER(l) 

L: ENDERl  =.TRUE. 

DELAY  =  L: ENDERl 
L: ENDERl  =  .FALSE. 

RETURN 

END 

C 

SUBROUTINE  PREPDCA 

+ 

COMMON  /QQDCP/DCASF (0:0) 

DCASF(O)  =  1.0/QDCASR(0)/D1MAX 

RETURN 

END 


m. 
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*PSP=1,0,BRR«ALL 

•TITLE 

EFFECTS  OF  FIRST  ORDER  ANTI-ALIASING  FILTER  ON  SA 

•INPUT 

PROGRAM 

•INTERRUPT  DECLARATIONS' 

•INTDEF(0,1,1) ' 

•INTDEF(1,1,0) ' 

INITIAL 

'eBETA(BETA) * 

MAXT  =  PERIOD/BETA 
LOGPER  =  Cl NT  *  BETA 
FRAD=FREQ^2^PI 
POLE=l/FILTER 
VARIABLE  TIME  =  0 


CONSTANT  BETA=1  ,  RUNTIM  =  3.0  ,  PERIOD=  .015625 


CONSTANT  CINT=. 015625,  FREQ=32,  PI=3. 14159 
CONSTANT  FILTER=38.4,K=1 
' §PARAMETER  BETA, RUNTIM, FREQ , FILTER, K ' 
•gMAXVAL  FREQ=  64,  FILTER=200,K=1000 ' 
'0MINVAL  FREQ=.01,  FILTER* . 01 , K=1 ' 

'0MAXVAL  BETA  =100,  RUNTIM  =50,  TIME=50' 
•eMINVAL  BETA=.l  ,  RUNTIM  =0  ,  TIME=0  ' 


NSTEPS  NSTEP=1 
D1  =  0 

END  $'OF  INITIAL' 

DYNAMIC 

'INTERRUPT  RATE  ERROR  DECLARATIONS' 
LOGICAL  ENDERl,  RATBRl,  ERRORl 
ENDERl  =  .FALSE. 

ERRORl  =  RATER 1 
DERIVATIVE 

CONSTANT  D1MAX=1 
CONSTANT  A1MAX=1 
CONSTANT  CMAX=1 
'0SCALE  D1=D1MAX' 

'0SCALE  A1=A1MAX' 

'0SCALE  C=CMAX' 

D1=A1 

'^PARALLEL' 

R=HARM(0,FRAD,0) 

A1=REALPL ( POLE, R, 0 ) 

C=D1 

TERMT(TIME  .GT.  RUNTIM) 
'eRECORD(REC01,,,,,,,,R,C,Al,Dl,,,, ) ' 
'DEFINE  INTERRUPT  CONTROL' 

GPIO  =  CLOCK (PERIOD) 

GPU  =  CLOCK  (LOGPER) 

'@INTRRT  1  =  GPIO' 

•eiNTRRT  2  =  GPU' 

RATER 1 =R ATERR (GPIO, ENDER 1 ) 

'0END  PARALLEL' 

END  $'OF  DERIVATIVE' 


5 


n 


END  $'0F  DYNAMIC 
TERMINAL 

END  $'0F  TERMINAL' 

END  $ ' OF  PROGRAM ' 

•TRANSLATE 

DCA(1)=D1 

PADC(1)=A1 

•OUTPUT 

•END 

SUBROUTINE  PREPl 

+ 

INCLUDE  El.SAFIL 
A1=QRPADC(0)^S:A1 
RETURN 
END 
C 

SUBROUTINE  POSTl 

+ 

INCLUDE  El.SAFIL 

COMMON  /QQDCP/DCASF(0:0) 

LOGICAL  DELAY 

CALL  QWDCAR(0^D1^DCASF(0) ) 

IF  (LrRATERl)  CALL  ZZRTERd) 
L:ENDER1  =.TRUE. 

DELAY  =  L:ENDER1 
L:ENDERI  =  .FALSE. 

RETURN 

END 

C 

SUBROUTINE  PREPDCA 
4- 

COMMON  /QQDCP/DCASF(0:0) 
DCASF(O)  =  1.0/QDCASR(0)/D1MAX 
RETURN 
END 
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*PSP=1,0,ERR»ALL 

•TITLE 

HARDWARE  SIMULATION  WITH  NOISE  AND  OUTPUT  FILTER 

•INPUT 

PROGRAM 

•INTERRUPT  DECLARATIONS' 

*INTDEF(0,1,1) ' 

'INTDEF(1,1,0) ' 

INITIAL 

'eBETA(BETA) ' 

MAXT  =  PERIOD/ BETA 
LOGPER  =  Cl NT  »  BETA 
FRAD=FREQ^2^PI 
POLE=l/FILTER 
VARIABLE  TIME  =  0 
•^EXTERN  NOISE' 

CONSTANT  BETA=1  ,  RUNTIM  =  3.0  ,  PERIOD=  .015625 
CONSTANT  CINT=. 015625,  FREQ=32,  PI=3. 14159 
CONSTANT  FILTER=38.4,K=1,  MNOISE=1.0 
•ePARAMETER  BETA, RUNTIM, FREQ, F I LTER, K ' 

'§MAXVAL  FREQ=  64,  FILTER=200,K=100 ' 

•eMINVAL  FREQ=.01,  FILTER= . 01 , K=1 ' 

•eMAXVAL  BETA  =100,  RUNTIM  =50,  TIME=50' 

•eSCALE  NOISE=MNOISE' 

•eMINVAL  BETA=.l  ,  RUNTIM  =0  ,  TIME=0  ' 

NSTEPS  NSTEP=1 
D1  =  0 

END  $'OF  INITIAL' 

DYNAMIC 

•INTERRUPT  RATE  ERROR  DECLARATIONS' 

LOGICAL  ENDERl,  RATERl,  ERRORl 
ENDERl  =  .FALSE. 

ERRORl  =  RATERl 
DERIVATIVE 

CONSTANT  D1MAX=1.2 
CONSTANT  FMAX=1.2 
CONSTANT  CMAX=1.2 
'@SCALE  D1=D1MAX' 

•esCALE  FILOUT=FMAX' 

'§SCALE  C=CMAX' 

PROCEDURAL (D1=A1) 

D1=A1 

END 

'8PARALLEL' 

R=HARM(0,FRAD,0) 

RN  =  R  -t-  NOISE 
A1=RN 

FILOUT=REALPL{POLE,D1,0) 

C=FILOUT^K 

TERMT(TIME  .GT.  RUNTIM) 

•8RECORD(REC01, ,,,,,, ,R,C,Al,Dl,NOISE,,, ) ' 
•DEFINE  INTERRUPT  CONTROL' 

GPIO  =  CLOCK (PERIOD) 


GPU  CLOCK (LOGPER) 
•eiNTRRT  1  =  GPIO' 

•gINTRRT  2  =  GPU' 
RATERlsRATERR (GPIO, ENDERl ) 
'8END  PARALLEL* 

END  $'OF  DERIVATIVE' 

END  $'OF  DYNAMIC 
TERMINAL 

END  $'OF  TERMINAL* 

END  $*OF  PROGRAM* 

•TRANSLATE 

DCA(1)=D1 
PADC(1)=A1 
CONNECT  AIN31=NOISE 
•OUTPUT 
•END 

SUBROUTINE  PREPl 

+ 

INCLUDE  El. TESTS 
A1=QRPADC(0)^S:A1 
RETURN 
END 
C 

SUBROUTINE  POSTl 

+ 

INCLUDE  El. TESTS 

COMMON  /QQDCP/DCASF ( 0:0) 

LOGICAL  DELAY 

CALL  QWDCAR(0,D1^DCASF(0) ) 

IF  (L:RATER1)  CALL  ZZRTER(l) 
L: ENDERl  =.TRUE. 

DELAY  =  L  .’ENDERl 
L: ENDERl  =  .FALSE. 

RETURN 

END 

C 

SUBROUTINE  PREPDCA 


COMMON  /QQDCP/DCASF(0:0) 
DCASF(O)  =  1.0/QDCASR(0)/D1MAX 
RETURN 
END 


•PSP^l^O^ERRi^ALL 

•TITLE 

THESIS  Sl.SATL- INPUT  SELECTOR/MONITOR 

•INPUT 

PROGRAM 

'THIS  PROGRAM  SIMULATES  A  DUAL  REDUNDANT,' 

'  ASYNCHRONOUS  SYSTEM  BY  FINDING  THE  DIFFERENCE  OF 
'  THE  INPUT  AND  THE  INPUT  DELAYED  ONE  SAMPLE.  ' 
'SINCE  THE  AFTI/F-16  RUNS  AT  256  SAMPLES/SECOND,' 

'  THE  SAMPLE  TIME  IS  .00390625  SEC.  THIS  SIMULATES 
'  THE  GREATEST  TIME  SKEW  BETWEEN  COMPUTERS' 

'OF  0.00390625  SECONDS' 

'INTERRUPT  DECLARATIONS' 

'INTDEF(0,1,1) ' 

'INTDEF(1,1,0) ' 

INITIAL 

'eBETA(BETA) ' 

MAXT  =  PERIOD/BETA 
LOGPER  =  Cl NT  *  BETA 
FRAD=FREQ*2*PI 
VARIABLE  TIME  =  0 

CONSTANT  BETA=1,RUNTIM  *  3.0,PERIOD=  .00390625 
CONSTANT  CINT=. 00390625,  FRAD=131 . 9 , FREQ=21 
CONSTANT  M=3,  PI=3. 14159 
'^PARAMETER  BETA, RUNTIM, FREQ ' 

'gMAXVAL  FREQ=256' 

'eMINVAL  FREQ=.01' 

•eMAXVAL  BETA  =100,  RUNTIM  =50,  TIME=50’ 
'@MINVAL  BETA=.l  ,  RUNTIM  =0  ,  TIME=0  ' 

NSTEPS  NSTEP=1 
D1  =  0 
D2=0 
DIF  =0 

END  $'OF  INITIAL' 

DYNAMIC 

•INTERRUPT  RATE  ERROR  DECLARATIONS' 

LOGICAL  ENDERl,  RATERl,  ERRORl 
ENDERl  =  .FALSE. 

ERRORl  =  RATERl 
DERIVATIVE 

CONSTANT  D1MAX=3 . 1, D2MAX=3 . 1 , DIFMAX=3 . 1 
CONSTANT  A1MAX=3.1 
•eSCALE  D1=D1MAX' 

'§SCALE  D2=D2MAX' 

•esCALE  DIF=DIFMAX' 

'@SCALE  A1=A1NAX' 

PROCEDURAL (D1,D2, DIF  =A1 ) 

'8IMPL(D1) • 

DIF=D2-D1 
D2  =  D1 
D1=A1 

'9END  IMPL' 

END 
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'^PARALLEL' 

R=HARM(0,FRAD,0) 

A1=M*R 

C=D1 

TERMT(TIME  .GT.  RUNTIM) 
'@RECORD(REC01,,,,,,,,R,DIF,D2,,,,, ) ' 
'DEFINE  INTERRUPT  CONTROL' 

GPIO  =  CLOCK (PERIOD) 

GPU  =  CLOCK  (LOGPER) 

'eiNTRRT  1  =  GPIO' 

'eiNTRRT  2  =  GPU' 

RATER1=RATERR (GPIO, ENDERl ) 

'@END  PARALLEL' 

END  $'OF  DERIVATIVE' 

END  $'OF  DYNAMIC 
TERMINAL 

END  $'OF  TERMINAL' 

END  $'OF  PROGRAM' 

•TRANSLATE 

DCA(1)=D1 
PADC ( 1 ) =A1 
•OUTPUT 
•END 

SUBROUTINE  PREPl 

+ 

INCLUDE  El.SATL 
A1=QRPADC(0)^S:A1 
RETURN 
END 
C 

SUBROUTINE  POSTl 

+ 

INCLUDE  El.SATL 
COMMON  /QQDCP/DCASF(0:2) 

LOGICAL  DELAY 

CALL  QWDCAR(0,D1^DCASF(0) ) 

CALL  QWDCAR(1,D2^DCASF(1) ) 

CALL  QWDCAR(2,DIF  •DCASF(2)) 

IF  (L:RATER1)  CALL  ZZRTER(l) 

L: ENDERl  =.TRUE. 

DELAY  =  L: ENDERl 
L: ENDERl  =  .FALSE. 

RETURN 

END 

C 

SUBROUTINE  PREPDCA 

+ 

COMMON  /QQDCP/DCASF (0:0) 

DCASF(O)  =  1.0/QDCASR(0)/D1MAX 
DCASF(l)  =  1.0/QDCASR(1)/D2MAX 
DCASF(2)  =  1.0/QDCASR(2)/DIFMAX 
RETURN 
END 
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*PSP=1,0,ERR=ALL 

‘TITLE 

THESIS  Sl.RLSIM  RATE  LIMITER 

‘INPUT 

PROGRAM 

'THIS  PROGRAM  SIMULATES  1:  THE  RATE  LIMITER' 

•  AFTER  THE  OUTPUT  S/M  WITH  NO  INTER-CHANNEL  ' 
•DIFFERENCE  BETWEEN  COMPUTERS  (RATE!),' 

•2:  THE  RATE  LIMITER  AFTER  THE  OUTPUT  S/M  WITH  ' 
'MAXIMUM  INTER-CHANNEL  DIFFERENCE  (RATE2)' 

'AND  3:  RATE  LIMITERS  PRIOR  TO  THE  OUTPUT  S/M  WITH' 
•MAXIMUM  INTER-CHANNEL  DIFFERENCE  (RATE3)' 

•THE  SAMPLE  PERIOD  IS  .015625  SEC,  WHICH  SIMULATES' 
'THE  AFTI  COMPUTATIONAL  FRAME  RATE' 

I  f 

'INTERRUPT  DECLARATIONS' 

' INTDEF(0,1,1) ' 

'INTDEF(1,1,0) ' 

INITIAL 

'eBETA(BETA) ' 

MAXT  =  PERIOD/BETA 
LOGPER  =  CINT  ‘  BETA 
FRAD=FREQ*2*PI 
VARIABLE  TIME  =  0 

CONSTANT  BETA=1,  RUNTIM=3.0,  PERIOD= . 015625 
CONSTANT  TL=.45,  PI=3. 14159 

CONSTANT  CINT=. 015625,  FRAD  =  131.9,  FRE0*21 
'ePARAMETER  BETA, RUNTIM, FREQ, TL ' 

'eMAXVAL  FREQ=256,  TL=1' 

•eMINVAL  FREQ=.01,  TL=0' 

'eMAXVAL  BETA  =100,  RUNTIM  =50,  TIME=50' 
'eMINVAL  BETA=.l  ,  RUNTIM  =0  ,  TIME=0  ' 

NSTEPS  NSTEP=1 
D1=0 
N=0 

END  $'OF  INITIAL' 

DYNAMIC 

'INTERRUPT  RATE  ERROR  DECLARATIONS' 

LOGICAL  ENDERl,  RATER 1,  ERROR 1 
ENDERl  =  .FALSE. 

ERRORl  =  RATERl 
DERIVATIVE 

CONSTANT  D1MAX=1,  R1MAX=1,  R2MAX=1,  R3MAX=1 
CONSTANT  A1MAX=1.0 
'eSCALE  D1=D1MAX' 

•eSCALE  RATE1=R1MAX’ 

'eSCALE  RATE2=R2MAX' 

'eSCALE  RATE3=R3MAX' 

'eSCALE  A1=A1MAX' 

PR0CEDURAL(RL1,RL2,RL3  =  D1,TL,N) 

D1=A1 

CALL  RL{RL1,RL2,RL3,D1,TL,N) 
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RATE1=RL1 

RATE2=RL2 

RATE3=RL3 

END 

'ePARALLEL' 

R=HARM(0,FRAD,0) 

A1=R 

C=D1 

TERMT(TIME  .GT.  RUNTIM) 

'eRECORD(REC01, , , , ,, , ,R,C,RATEl,RATE2,RATE3, 
•DEFINE  INTERRUPT  CONTROL’ 

GPIO  =  CLOCK (PERIOD) 

GPU  =  CLOCK  (LOGPER) 

•ilNTRRT  1  =  GPIO* 

•eiNTRRT  2  =  GPU' 
RATER1=RATERR(GPI0,ENDER1 ) 

’SEND  PARALLEL' 

END  $'OF  DERIVATIVE* 

END  $'OF  DYNAMIC 
TERMINAL 

END  $'OF  TERMINAL* 

END  $*OF  PROGRAM* 

•TRANSLATE 

DCA(1)=D1 

PADC(1)=A1 

•OUTPUT 

•END 

SUBROUTINE  PREPl 

+ 

INCLUDE  El.RLSIM 
A1=QRPADC(0)*S:A1 
RETURN 
END 
C 

SUBROUTINE  POSTl 

+ 

INCLUDE  El.RLSIM 
COMMON  /QQDCP/DCASF(0:2) 

LOGICAL  DELAY 

CALL  QWDCAR(0,D1*DCASF(0) ) 

CALL  QWDCAR(1,RATE1^DCASF(1) ) 

CALL  QWDCAR( 2,RATE2  *DCASF(2)) 

CALL  QWDCAR(3,RATE3  *DCASF(3)) 

IF  (L:RATER1)  CALL  ZZRTER(l) 

L:ENDER1  =.TRUE. 

DELAY  =  L:ENDER1 
L:ENDER1  =  .FALSE. 

RETURN 

END 

C 

SUBROUTINE  PREPDCA 


+ 


COMMON  /QQDCP/DCASF(0:0) 


DCASF(O)  =  1.0/QDCASR(0)/D1MAX 
DCASF(l)  =  1.0/QDCASR(1)/R1MAX 
DCASF(2)  =  1.0/QDCASR(2)/R2MAX 
DCASFO)  =  1.0/QDCASR(3)/R3MAX 
RETURN 
END 

SUBROUTINE  RL(R1,R2,R3, Zl,TLEV,C) 

IF  (C  .EQ.  0)  THEN 
R  =  0 
R1  =  0 
R2  =  0 
Z2  =  Z1 
END  IF 
C=C+1 

DIF1=R1-Z1 

IF  (DIFl  .LE.  -TLEV)  THEN 
R1=R1+TLEV 

ELSE  IF  (DIFl  .GE.  TLEV)  THEN 
R1=R1-TLEV 
ELSE 

R1  =  Z1 
END  IF 

AVE2=(Z2+Zl)/2 

DIF2=R2-AVE2 

IF  (DIF2  .LE.  -TLEV)  THEN 
R2=R2+TLEV 

ELSE  IF  (DIF2  .GE.  TLEV)  THEN 
R2=R2-TLEV 
ELSE 

R2*AVE2 
END  IF 
DIF  =  R-Z2 

IF  (DIF  .LE.  -TLEV)  THEN 
R=R+TLEV 

ELSE  IF  (DIF  .GE.  TLEV)  THEN 
R=R-TLEV 
ELSE 
R  =  22 
END  IF 
R3=(Rl+R)/2 
Z2  =  Z1 
RETURN 


Appendix  C;  Resources  Required 

The  following  is  a  list  of  resources  used  to  simulate 
nd  study  subharmonic  aliasing. 

a)  Cyber  Computer  (Total  CAD  Package) 

b)  ICECAP  CAD  Package 

c)  SIMSTAR  Hybrid  Computer 

d)  Strip  Chart  Recorder 
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result  of  impostor  frequencies  (much  like  aliasing)  being 
introduced  into  the  output  signal  by  the  sampling  process. 

Well  defined  subharmonic  aliases  occur  due  to:  impostor 
frequencies  and  a  phenomena  known  as  apparent  low  frequency 
surge,  which  occurs  when  the  input  frequency  is  nearly  an 
integer  multiple  (>1)  of  the  impostor  frequency.  “7^ 

Simulations  demonstrate  the  effectiveness  of  input  ^ 

analog  low  pass  filters  for  attenuating  signals  in  this 
frequency  range.  Unfortunately,  according  to  simulations, 
interchannel  differences  may  exceed  Input  or  Output  selector 
monitor  threshold  levels  for  frequencies  in  the  subharmonic 
range.  Simulations  also  show  that  it  is  the  combination  of 
well  defined  subharmonic  aliases,  the  high  frequencies 
(relative  to  the  sample  rate)  of  the  input  signal,  and  the 
nonlinear  characteristics  of  the  software  rate  limiter  which 
form  the  low  frequency  output  from  the  rate  limiter. 

In  conclusion,  since  subharmonic  aliasing  occurs  for 
frequencies  which  are  high  relative  to  the  sample  rate, 
careful  selection  of  the  sample  rate  (and  computational 
frame  rate)  should  reduce  this  undesireable  effect. 

Additionally,  since  the  phenomenon  is  actually  the 
combination  of  two  high  frequencies,  the  use  of  analog 
lottfpass  filters  at  the  input  does  attenuate  the  subharmonic 
aliasing  effect. 
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